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(54) Computation of coefficients for a decision feedback equaliser with variable delay 



(57) Decision Feedback Equalizer (DFE) coeffi- 
cients are determined from a channel estimate by cast- 
ing the DFE coefficient problem as a standard recursive 
least squares (RLS) problem and solving the RLS prob- 
lem. In one embodiment, a fast recursive method, e.g., 
fast transversal filter (FTP) technique, is used to com- 
pute the Kalman gain of the RLS problem, which is then 
directly used to compute MIMO Feed FonA/ard Equalizer 



(FFE) coefficients. The FBE coefficients are computed 
by convolving the FFE coefficients with the channel im- 
pulse response. Complexity of a conventional FTF al- 
gorithm may be reduced to one third of its original com- 
plexity by selecting a DFE delay to force the FTF algo- 
rithm to use a lower triangular matrix. The length of the 
DFE may be selected to minimize the tap energy in the 
FBE coefficients or to ensure that the tap energy in the 
FBE coefficients meets a threshold. 
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Description 

CROSS REFERENCES TO RELATED APPLICATIONS 

5 [0001] This application is a continuation-in-part of U.S. Regular Utility Application Serial No. 10/154,947, filed May 
24, 2002 which clainns priority to U.S. Provisional Application Serial No. 60/339,253, filed October 26, 2001 and which 
is a continuation-in-part of U.S. Regular Utility /Vpplication Serial No. 10/044,013, filed October 26, 2001 , which claims 
priority to U.S. Provisional Application Serial No. 60/322,994, filed September 1 8. 2001 . This application clainns priority 
to each of these applications and incorporates herein each of these references in its entirety for all purposes. 
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1. FIELD OF THE INVENTION 

[0002] This invention relates generally to digital communications; and more particularly to decision feedback based 
equalizers that are employed in digital communication systems. 

75 

2. BACKGROUND OF THE INVENTION 

[0003] The structure and operation of communication systems is generally known. Many communication systems 
carry data, e.g., voice, audio, video, file, or other digital data that is sent from a transmitter to a receiver. On the trans- 
it? mitter side, data is first formed into packets. This data may be raw data or encoded data that represents the raw data. 
Each of these packets also typically includes a header, a known training sequence, and a tail. These packets are then 
modulated into symbols and the symbols are transmitted by the receiver and intended for receipt by the receiver. The 
receiver then receives the symbols and attempt to extract the data from the packets that are canied by the symtK3ls. 
[0004] A "channel" carries the symbols from the transmitter to the receiver. The channel is serviced by a wired, 
25 wireless, optical, or another media, depending upon the communication system type. In many communication systems, 
such as terrestrial based wireless communication systems, satellite based communication systems, cable based com- 
munication systems, etc., the channel distorts the transmitted symbols, from the perspective of the receiver, causing 
interference between a subject symtjol and a plurality of symbols surrounding the subject symbol. This type of distortion 
is referred to as "inter-symbol-interference" andMs, generally speaking, the time-dispersed receipt of multiple copies 
30 the symbols caused by multipath. The channel also introduces noise into the symbols prior to their receipt. Each of 
these concepts is well known. 

[0005] Equalizers are now generally employed in an attempt to remove channel effects from a received symbol 
stream. Thus, equalizers are essential building blocks of modern receivers, especially in broadband applications where 
inter-symbol-interference is a critical problem. In a typical equalizer, the channel between the transmitter and the re- 

35 ceiver is first estimated based upon the training sequence contained in one or more preambles. Then optimal equalizer 
coefficients (also referred to as taps and/or tap coefficients for the equalizer) are estimated based upon the channel 
estimate. The optimal equalizer coefficients are then used by the equalizer in extracting the data from the packet. The 
optimal equalizer coefficients may also be computed after extraction of the data from the equalized data stream based 
upon blind channel estimates. Equalizer coefficient generation should be repeated as often as possible, especially in 

40 fast varying channel cases, to generate new equalizer coefficients. The received data stream is usually buffered during 
the period that is required for channel estimation and equalizer coefficient computations. Thus, the preamble (and also 
actual data) contained in a packet may be used to generate the channel estimate and optimal equalizer coefficients 
that are employed by the equalizer to extract the data from the packet. 

[0006] As symbol rates increase and modulation schemes become more complex, equalizers have increasingly 
45 greater importance. A critical factor in increasing the effectiveness of these equalizers is the complexity of optimal 
equalizer coefficient computation. A reduction in this complexity: (1) reduces the memory size required to buffer the 
received symbol stream sequence during the period required for coefficient computations: (2) allows more frequent 
uploading of new coefficients thus enabling the equalizer to track fast channel variations; and (3) simplifies the hardware 
and, resultantly, the die area required for coefficient computation. 
50 [0007] FIG. 1 is a block diagram illustrating a discrete time symbol-spaced Decision Feedback Equalizer (DFE) based 
channel equalization model 100. The channel equalization model 100 includes a channel 102, a Feed Fonward Equal- 
izer (FFE) 104, a Decision block 106, and a Feed Back Equalizer (FBE) 108. An input sequence x(n) is complex, 
independent and identically distributed with unit power. Additive noise v(n) is white Gaussian with power spectral density 
o^. Furthermore, the decisionsjt (rhb) are assumed to be correct, and hence equal to x(n - 5). This assumption makes 
55 the design of the FBE 108 and FFE 104 easier, but at the expense of introducing error propagation due to possibly 
wrong decisions. The FFE 1 04 function G(z) has length L The channel (impulse) response vector of the channel h is 
given in Equation (1) as: 
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Equation (1) 



5 The number of coefficients (taps) M of the FBE 1 08 function B(z) is assumed greater or equal to the channel memory, 
i.e., M>N'1. These modeling assumptions are feasible in practice. 

[0008] In estimating FFE 104 and FBE 108 equalizer coefficients, the goal is to minimize the mean square error 
quantity of Equation (2). 
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C = E|x(n.5)-x(n-5)r. 



Equation (2) 



where x(/>8)is the delayed input signal estimate prior to the Dejjision block 106. By collecting the coefficients of both 
G{z) and 6(z) into vectors, we can express the received signal x(n - 6) in Equation (3) as: 



x^^yng-x^b Equation (3) 

A channel output model defining y„ may be expressed by: 

y„=x„H+v„ Equation (4) 

where H is the (A/ + L -1 ) x I convolution matrix corresponding to the channel response and expressed as: 
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In this model, is the 1 x (A/ + L-^ ) input vector, 



Equation (5) 
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y„ is the 1 X L input regression vector to the FFE 104, 



Equation (6) 



Equation (7) 



X n is the 1 X M input regression vector to the (strictly causal) FBE 108. 
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x„L[xin-S-\) x{n-S-2) jc(/i-^-M)] Equatioii(8) 

5. 

and v„ is the 1 x /. vector noise process. 

[0009] The current efficient methods for computing the optimal filter coefficients of a decision feedback equalizer, 
which optimizes (2), are based on the well-known Choiesky decomposition method (from a finite-dimension problem 
fomnulation). Two published papers: (1) N. Al-Dhahir and J.M. Cioffi, "MMSE Decision-Feedback Equalizers: Finite- 
^0 Length Results," IEEE Trans, on Infomnation Theory, vol. 41 , no. 4. pp. 961 -973, July 1 995; and (2) N. Al-Dhahir and 
J.M. Cioffi, 'Fast Computation of Channel-Estimate Based Equalizers in Packet Data Transmission," IEEE Trans, on 
Signal Processing, vol. 43, no. 11 , pp. 2462-2473, Nov. 1 995 provide one procedure for computing optimal DFE settings. 
These equations are referred to hereinafter as "Al-Dhahir's equations." 

[001 0] Generally speaking, Al-Dhahir's equations, as well as other existent techniques rely on the use of the gener- 
J5 alized Schur algorithm for fast Choiesky decomposition of the matrices Involved in both the FBE and FFE optimal 
coefficient setting computations. However, the overall procedures for calculation of the DFE (FBE and FFE) coefficients 
have the folbwing problems: 

1 . These procedures require the use of nonstructured recursive equations. These equations are difficult to imple- 

20 ment from the perspective of integrated circuit design. In particular, the recursive equations used in Al-Dhahir's | 

DFE tap (coefficient) computer often requires the use of a DSP processor. As Is generally known, the use of a 
DSP processor in a real-time communication system application severely constricts system throughput. 

2. These procedures, in particular Al-Dhahir*s equations, require a complex tap computer In other words, they 
require the use of a relatively large number of complex multiplies. 

25 3. These prior art DFE coefficient computation techniques can be used only for equalizers that use a fixed equalizer 

delay (5), which is set to its maximum value 1-1 , where L is the length of the FFE 104. The prior art techniques 
cannot use a different delay 5 of the equalizer. 

[001 1 ] Additional difficulties relate to the operation of FFEs. and the computation of FFE coefficients where the chan- 
30 nel being equalized has a long pre-echo (precursor components). The problem arises from the modeling restriction 
that the total channel/equalizer delay be less than the FFE length. White this restriction enables the derivation of fast 
computation of the FBE taps, such as with Al-Dhahir's equations, this restriction results in FBE taps greater than one, 
which causes the DFE to be unstable. 

[0012] Thus, there is a need in the art for an equalizer coefficient computation methodology and system in which 
35 such methodology may be implemented. 

SUMMARY OF THE INVENTION 

[0013] In order to overcome the above-cited shortcomings, among others of the prior art, the method and system of 

40 the present invention computes optimal Decision Feedback Equalizer (DFE) coefficients {goppb^pi) from a channel 
estimate h. As contrasted to the prior art, the methodology of the present invention relies on iterating simple structured 
recursions to yiekl the DFE coefficients. Resultantly, the operation according to the present invention produces DFE 
coefficients faster and with less computational requirements than the prior art operations, thus providing significant 
advantages over prior techniques. 

45 [001 4] According to the present invention, a channel impulse response h is first estimated based upon either a known 
training sequence or an unknown sequence. A solution to the DFE coefficient computation problem is then cast as a 
standard recursive least squares (RLS) problem. More specifically, the solution for Feed Forward Equalizer (FFE) 
coefficients g^pf (of the DFE) based upon the channel Impulse response h Is formulated as the Kalman gain solution 
to the RLS problem. A fast recursive method for computing the Kalman gain is then directly used to compute g^pf. 

50 [001 5] In one embodiment of the present invention, a fast transversal filter (FTF) technique is employed to compute 
g^pf. In this embodiment, the complexity of a conventional FTF algorithm is reduced to one third of its original complexity 
by choosing the length of Feed Back Equalizer (FBE) coefficients b^pf (of the DFE) to force the FTF algorithm to use 
a lower triangular matrix. This technique significantly reduces the needed recursions and computations, as well as 
avoiding finite precision problems of a conventional FTF solution implemented in hardware. 

55 [0016] Finally, with the FFE coefficients g^pf determined, the FBE coefficients b^p^ are computed by convolving the 
FFE coefficients with the channel impulse response h. In perfomning this operation, a convolution matrix that 
characterizes the channel impulse response h is extended to a bigger circulant matrix. With the extended circulant 
matrix structure, the convolution of the FFE coefficients dopr with the channel Impulse response h may be performed 
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the convolution operations in a transformed domain. In one embodiment, the convolution is performed in the frequency 
domain, which can be computed efficiently using the Fast Fourier Transfomi (FFT). However, in other embodiments, 
other transfonmation domains are employed for the convolution operations, e.g., Discrete Cosine Transfonnation do- 
main and the Discrete Hadamard Transfomn domain, among others. 

5 [0017] According to one aspect of the present invention, solving the least squares problem to yield Feed Fonvard 
Equalizer (FFE) coefficients of the DFE coefficients includes: (1) selecting a delay of the DFE; (2) when the delay is 
less than or equal to a number of taps of the FFE minus one, performing a first set of operations; and (3) when the 
delay is greater than the number of taps of the FFE minus one, performing a second set of operations. Generally, this 
aspect of the present invention Irees' the restriction of the channel/equalizer delay, enables a fast DFE tap computation 

^0 algorithm, and also results in a FBE coefficient solution with a long delay (longer than the FFE length). This technique 
is perfonned at the lowest possible increase in complexity in the DFE tap computation algorithm. With a delay greater 
than the FFE length, it turns out that the computation iterates forward prediction recursions as well as backward pre- 
diction recursions. The algorithm is recursive with D iterations (i = 1 to D). For i < FFE length (L), only forward prediction 
recursions are iterated. For i = L to D, both forward and backward recursions are iterated. For the first L iterations, 3L 

15 multiplies/iteration are required. For the last (D-L) iterations, the complexity is 5L multiplies/iteration. Thus the total 
complexity is [SL^ + 5L (D-L)) multiplies. 

[0018] With the constraint of fixed FBE delay removed, however, it is now required to select a FBE delay. Thus, 
another aspect of the present invention focuses on determining the FBE delay. According to one operation of the 
present Invention determining the delay of the DFE includes determining a delay of the channel based upon the channel 
20 response of the channel and adding the delay of the channel to the number of taps of the FFE minus one to produce 
the delay of the DFE. In this operation, the delay of the channel may be selected as the length of the precursor of the 
channel response. 

[0019] Another technique for detenmining the delay of the DFE according to the present invention considers the tap 
energy in the FBE coefficients. One embodiment of this technique includes: (1) calculating a plurality of sets of FFE 

25 coefficients and FBE coefficients, wherein each of the plurality of sets of FFE coefficients and FBE coefficients corre- 
sponds to a respective DFE delay and (2) selecting one of the plurality of sets of FFE coefficients and FBE coefficients 
based upon the tap energy of the respective FBE coefficients. With this embodiment, in selecting one of the plurality 
of sets of sets of coefficients includes selecting a set of FFE coefficients and FBE coefficients that meets an FBE 
coefficient tap energy threshold. In an alternate operation, selecting one of the plurality of sets of sets of coefficients 

30 includes selecting a set of FFE coefficients and FBE coefficients that provides a minimum FBE coefficient tap energy. 
[0020] The method of the present invention is exact and is much faster than any prior technique that may have been 
employed to compute DFE coefficients for the same channel and DFE filter lengths. Because the method of the present 
invention relies upon simple structured recursion, the computation period and/or the hardware required for the DFE 
coefficient computations is significantly reduced as compared to prior art implementations. The reduction in computation 

35 complexity and resultant increase in computation speed enables a DFE that operates according to the present invention 
to track fast channel variations. 

[0021] The method of the preset invention may be efficiently used in many communication systems that require the 
use of equalizers such as mobile wireless receivers, fixed wireless receivers, cable modem receivers, MDTV receivers, 
etc. to increase the perfomnance of such communication systems. Other features and advantages of the present in- 
^0 vention will become apparent from the following detailed description of the invention made with reference to the ac- 
companying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

^5 [0022] These and other features, aspects and advantages of the present invention will be more fully understood 
when considered with respect to the following detailed description, appended claims and accompanying drawings 
wherein: 

FIG. 1 is a block diagram illustrating a discrete time symbol-spaced Decision Feedback Equalizer (DFE) channel 
50 equalization model 100; 

FIG. 2 is a logic diagram generally illustrating operation according to the present invention in determining DFE 
coefficients and applying such coefficients to a DFE; 

FIG. 3 is a logic diagram illustrating operations according to the present invention employed to detemriine Feed 
Forward Equalizer (FFE) coefficients for the DFE; 
55 FIG. 4 is a logic diagram illustrating operations according to the present invention employed to determine Feed 

Back Equalizer (FBE) coefficients for the DFE; 

FIG. 5 is a block diagram illustrating a discrete time fractionally-spaced DFE that operates according to the present 
invention; 
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FIG. 6 is a block diagram illustrating a multi-channel equivalent of the discrete time fractionally-spaced DFE of 
FIG. 5; 

FIG. 7 is a block diagram illustrating a transceiver constructed according to the present invention; 
FIG. 8 is a block diagram illustrating a Multi-Input-Multi-Output (MIMO) digital communication system that operates 
5 according to the present invention to equalize a channel; 

FIG. 9 is a system diagram illustrating a wireless digital. communication system in which a MIMO receiver 904 
operates according to the present invention; 

FIG. 10 is a system diagram illustrating a wireless digital communication system that includes a plurality of trans- 
mitters and a MIMO receiver that operates according to the present invention; 
10 FIG. 1 1 is a system diagram illustrating a wired digital communication system that includes a plurality of transmitters 
and a MIMO receiver that operates according to the present invention; 

FIG. 12 is logic diagram illustrating operatbn according to another embodiment of the present invention in deter- 
mining DFE coefficients; 

FIG. 13 is a graph illustrating a hypothetical channel impulse response indicating how channel delay may be 
IS detemnined according to the present invention; 

FIG. 1 4 is a logic diagram illustrating operation according to a first embodiment of the present invention in selecting 
FBE coefficients; and 

FIG. 15 is a logic diagram illustrating operation according to a second embodiment of the present invention in 
selecting FBE coefficients. 

20 

DETAILED DESCRIPTION 

[0023] FIG. 2 is a logic diagram generally illustrating operation according to the present invention in detenmining 
Decision Feedback Equalizer (DFE) coefficients and in applying such coefficients to a DFE. The operations of the 

25 present invention are performed by a procesisor, such as a Digital Signal Processor (DSP), or other circuitry present 
within a receiver that detemnlnes DFE coefficients to be applied to a DFE, also resident in the receiver. The DFE 
operates upon samples of a received signal in an attempt to remove channel effects from the samples so that digital 
data may be extracted from the samples. The structure and operation of DFEs, one of which was illustrated in FIG. 1 , 
are generally known and will not be further described herein except as they relate to the present invention. 

30 [0024] A processor, tap computer, DSP, or other receiver device, to detennine initial DFE coefficients to be used in 
subsequent operations by the receiver, will first perform the operations of the present invention. Thus, during startup 
or reset, a channel corresponding to the receiver is estimated (step 202). According to one embodiment of the present 
invention, the channel is estimated based upon a known preamble sequence. However, in other embodiments, the 
channel could also be estimated based upon unknown received data. In either case, channel estimation operations 

35 are generally well known and are not described further herein except as it relates to the present invention. 

[0025] With the channel estimated, Feed Forward Equalizer (FFE) coefficients are determined based upon the chan- 
nel estimate (step 206). Then, Feed Back Equalizer (FBE) coefficients are determined based upon the FFE coefficients 
and the channel estimate (step 208). The manner in which the FFE and FBE coefficients are generated is step 206 
and step 208 will be described in detail herein with reference to FlGs. 3-7. 

40 [0026] With the FFE and FBE coefficients determined, they are applied to the DFE (step 208) and are used In equal- 
izing samples of a received signal to remove channel effects. These DFE coefficients are continually updated (step 
210) using a known technique. Periodically, upon the receipt of a next packet for example, upon an indication that a 
new determination is required, or upon another triggering event, the DFE coefficients are again determined (step 212). 
In this event, another channel estimate may be obtained (step 21 4). Then, the DFE coefficients are again determined 

45 according to the present Invention and applied to the DFE. The operations of FIG. 2 continue until the receiver is turned 
off, placed in a sleep mode, or othenwise inactivated. 

[0027] FIG. 3 is a logic diagram illustrating operations according to the present invention employed to determine 
Feed Fonvard Equalizer (FFE) coefficients for the DFE. In a first operation of FIG. 3, a DFE delay is selected (step 
302). In one embodiment, this delay is selected as the channel length. In such case, the DFE delay corresponds to 
so the length of the FFE. 

[0028] Next, the DFE solution is formulated into a least squares solution (step 304). By collecting the FFE coefficients 
g and the FBE coefficients b into a single vector w, the minimization of the Mean Square Error may be written as: 
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minE 



Equation (9) 



[0029] Now, denoting R^J the variance of the augmented input regression vector u, and cross variance Rux(fhSi* ^® 
well-known solution to this snnoothing problem is given by Equation (10) as: 



Equation (10) 



where 



i?„=E 



yn 



Equation (11) 



and 



Equation (12) 



[0030] Using the channel output model of Equation (4), and the fact that x(n) is individually identically distributed (i. 
i.d.), the following closed form expressions for 



are detemnined: 



R^ = Ey„y„ = aJ/+H*H 



Equation (13) 



tn: 



Equation (14) 



Equation (15) 



R^^,„.,,=H'Ex:x{n-S)^H%^, 1 Oj ^h' 



Equation (16) 



EP 1 365 554 A1 



= 0 Equation (17) 

where H is a submatrix of H as set forth in Equation (5)» 



H 

L^2J 



Equation (18) 



is defined as the (5+1) x L submatrix of H, consisting of the first (5 + 1 ) rows of H. Note that for the case of colored 
noise, the matrix o^/ should be replaced by the autocon^elation matrix of the noise Extending the derivation to the 
colored noise case is straightforward. 

[0031] Now, with the quantities defined above, Equation (10) becomes: 



Using the well known inverse formula of block matrices, w^p^ may be rewritten as: 



^-.^=[[0 J] + 1^ ^ + -//•//)- [/ Equation (20) 



and as: 



which may be written as: 



= j^^j(<7,V + //;//, Equation (21) 



[0032] Although each of the two matrices Hj and has a shift structure, the augmented matrix 



does not have a shift structure. 

[0033] In selecting the length of the FBE (DFE delay) to be Af ^ A^l , the quantities and H2 are such that: 
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H,\ [0 Hj 



Equation (23) 



and 



Equation (24) 



This implies that: 



15 



Equation (25) 



2^ is block diagonal. In this case, the expressions for and b^p, decouple into a simple form. Therefore, the optimal 
FFE and FBE coefficients are represented as: 



25 



Equation (26) 



30 



Equation (27) 



[0034] The above expressions are valid for all values of the DFE delay 6. In general, the optimal value for the DFE 
delay 5 is within the range H <5op^A/+L-2. In the special case of the choice 5 = L-1 , the matrices involved in the above 
expressions are given by 



35 
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• 0 




' h'{N-\)' 






hiO) ■ 


• 0 




h'{N~2) 




A(Ar-i) 


h{N-2) ■ 


. A(0). 




. A*(0) . 



Equation (28) 



and 
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0 h{N-\) h{\) 
0 0 h(2) 



0 



h{N-l) 



Equation (29) 



[0035] Note that this solution is equivalent to the solution of Al-Dhahir's equations, due to the uniqueness of vv^pj 
when Equation (2) is minimized. However, the expressions obtained above for Equations (26) and (27) provide alter- 
^ native methods to compute the FFE and FBE coefficients and in a simpler and more efficient way than the 
prior art techniques. 

[0036] in calculating p^pj, a coefficient matrix is first defined in Equation (30) as: 
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Ps = {(^ll + H^^H^y Equation (30) 

5 

SO that the optimal solution for the FFE coefficients g^^x given by 

^op*=P5^5- Equation (31) 

w 

[0037] The expression for Qqp^ con'esponds exactly to the definition of the Kalman gain vector, used to update the 
optimal weights in a certain regularized Recursive Least Squares (RLS) problem. More specifically, given a (n + 1) x 
L data matrix H„ and the corresponding coefficient matrix P^, the Kalman gain = Pfp can be time-updated according 
to the following recursions: 

f^n)=1+^Pn-7^» Equation (32) 

20 On^Pfi'M") . Equation (33) 

Pn ='Pn-r9ny^ i")9„^ Equation (34) 

where P., = a ^ / (the term P.^ is the initial condition for P„) and gQ = 0. The computation of the Kalman gain vector p^^, 
in the above solution relies on the propagation of the Riccati variable P„. This method of computation requires O(L^) 
operations per iteration. 

[0038] Well known fast RLS schemes avoid the propagation of P„ and compute the gain g„ in a more efficient way. 
In this case, the computational complexity required is only of 0(L) per iteration. This implies that the overall complexity 
^ needed to calculate the FFE coefficients is on the order of 0{L^) thus yielding an efficient method for computing the 
FFE coefficients. Therefore, fast RLS filters are used to determine g^pf. In such case, fast transversal computations 
are employed to determine the Kalman gain for the RLS solution (step 306). Here we also note that it is straightfonward 
to extend this method to use other fast RLS algorithms, e.g., array form fast RLS algorithms. 
[0039] Faster recursions may be obtained by selecting the FBE length. Such selection eliminates certain variables 
that remain constant and equal to zero during adaptations in the special solution for g^pf. Fast RLS recursion in its 
explicit form propagates p„ in an efficient manner, see, e.g.: [1] Ljung, M. Morf, and D. Falconer, "Fast calculation of 
gain matrices for recursive estimation schemes," InL J. Contr. vol. 27, pp. 1-19, Jan 1978); [2] G. Carayannis, D. 
Manolakis, and N. Kalouptsidis, "A fast sequential algorithm for least squares filtering and prediction," IEEE Trans, on 
Acoustic. Speech, Signal Proc, vol. ASSP-31, pp. 1394-1402, December 1983; and [3] J. Cioffi and T. Kailath, "Fast 
recursive-least-squares transversal filters for adaptive filtering," IEEE Trans, on AcousL, Speech Signal Processing, 
vol. ASSP-32, pp. 304-337, April 1984. 

[0040] Table 1 lists the fast recursions employed in one embodiment for computing the normalized Kalman gain. The 
additional index included in /C/^ „ and yf^(n) stems from the fact that these quantities admit an order-update relation, 
instead of time-update. 
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Initialization 

>»^i=M'J =*o =0 

r.(0)=i 



Forn = 0 to Srepeat operations (1) to nS): 



(1) 
(2) 

(3) 

(4) 
(5) 

(6) 

(7) 
(8) 
(9) 

(10) 

(II) 
(12) 
(13) 



a,(n-l) = hin)-h„_,v{_^ 
/(«-l) = yi(«-l)a(n-l) 
a*(n-l) 



CHn-\) = CHn-2)+a*in- !)/(« - 1) 

y,.,(«)=y,(>i-l) 

v(n) = (last entry of 

*i^=*w^-iO:^) + v(nK-i 
;5(«) = ^»(«-l)v(«) 

n(«) = - 



i-r..,(«)>5.(«K(n) 

fc(n) = 7'i('«)A(«) 

^*(n) = C*(n-l) + ^*(«)K«) 



Table 1: Fast Transversal Computation of the Kalman gain. 



[0041] The purpose of computing k„ in the well-l<nown FTP algorithm is to use it in the computation of the corre- 
sponding optimal least squared solution for the FFE coefficients. Because we are interested only in k„, the filtering part 
of the FTF algorithm is not necessary, and does not appear in the algorithm listing. 

[0042] The quantities [w', vJ'] are known as the least-squares solutions of the forward and backward prediction 
problems of order L, with corresponding residual errors {/(n), b{n)}. Now, because during the first L iterations, the 
desired signal for the backward prediction problem is equal to zero, the backward least squares solution kv** will be 
equal to zero. This means that all quantities associated with the backward prediction problems will remain nul& for the 
first L iterations. Since in our case the optimal solution is achieved exactly at the L-th iteration, we can simply rule out 
the computation of these quantities from Table 1 . These operations correspond to operations (7) and (9) through (13) 
of Table 1 . Note that operation (1 0) of Table 1 is also eliminated, since ^(n) = 0, which implies Yl (n) = Y{.+i (n). This 
means that we can replace operatkxi (6) of Table 1 with: 



C'(l7-1) 



Equation (35) 



which can be further simplified since 
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yjn) = ^U^Mj^B:^^^ Equation (36) 

Moreover, operation (8) of Table 1 becomes simpty 

^L,n^f<L^^,n.^O'^)' Equation (37) 

where (1 :L) denotes the first L entries of 

[0043] Table 2 illustrates a simplified fast recursion for computing the optimal FFE coefficients g^p^ 



Initialization 



For n-Oto 5 repeat operations fl) through (7): 
(2) /(«-!) =y{/t-l)a(«-l) 



(3)* 



(4) ^^(n-l) = C^(/i-2) + a*(«-l)/(n-l) 
(6)K«) = - 



(7) *..=Wia-i) 



Table 2: Fast Transversal Computation of the FFE coefficients. 



[0044] With the Kalman gain determined, the FFE coefficients are then detennined (step 308). In the recursions of 
Table 2, the FFE coefficients are determined by setting g^pf = ki^y(b) when the number of iterations, n, is equal to the 
DFE delay. 

[0045] Unlike the conventional weighted FTF algorithm, the above recursions do not face finite precision difficulties 
for the following reasons. First, by ruling out the equations associated with the backward prediction problem, we are 

automatically eliminating many of the recursive loops that are responsible for the finite precision difficulties of the full 
FTF algorithm. Second, these simplified fast recursions have a forgetting factor X = 1, which yields finite precision 
stability. Third, the simplified algorithm deals with a finite set of data (5 +1 iterations). This algorithm is then reset, which 
avoids the accumulation of finite precision errors. 

[0046] FIG. 4 is a logic diagram illustrating operations according to the present invention employed to determine 
Feed Back Equalizer (FBE) coefficients for the DFE. The FBE coefficients d^p^ are determined according to the matrix- 
vector product of Equation (27). The computation of the feedback filter coefficients simply amounts to convolution of 
the channel impulse response with g^f. The convolution operation that defines the optimal FBE coefficients in Equation 
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(27) can be computed directly with LM/2 multiplications. Alternatively, the operations may also be efficiently performed 
using well-known fast FFT convolution techniques. 

[0047] To illustrate such operations. Equation (29) is first rewritten in Equation (38) as: 

^.=1/^ 0][^^ ^!][^o"]' Equation(38) 

The above equality holds true regardless of the values of H3, H4. or H5. Now, If these values are chosen such that the 
matrix O, defined by: 



L^4 Hy 



Equation (39) 



is a square (O x O) circulant matrix, where Q is the smallest power-of-two integer larger than or equal Mf L. In this 
case, the matrix C is rewritten in Equation (40) as: 

C = PAF, Equation (40) 

where F is a (0 x O) FFT matrix and A is a diagonal matrix that contains the elements of the FFT of the first row of 
C. The solution for b^pf becomes: 

*^=l^M 0]F*AFpJJ"j. Equation (41) 

35 The complexity of this method is the complexity of obtaining the FFT of two vectors of Q elements of each, the inverse 
FFT of another O element vector, and Q complex multiples. Thus the overall complexity is 



0 + 30 log2 (O). Equation (42) 

For the case of a power-of-two channel estimate N, the complexity is 

2M + 6M logs (2M), Equation (43) 



[0048] Thus, referring again to FIG. 4, in determining the FBE coefficients, the convolutional matrix H is first deter- 
mined (step 402). As was previously described, the matrix H may be detemnined as part of the channel estimation. 
Then, the convolutional matrix H is extended to form the bigger circulant matrix C (step 404). The bigger circulant 
matrix C is then converted to the Frequency domain (step 406) as are the FFE coefficients (step 408). With both the 
bigger circulant matrix C and the FFE coefficients in the frequency domain, they are convolved in the frequency domain 
by simple matrix multiplication (step 410) to produce the FBE coefficients b^p^ Finally, the resultant FBE coefficients 
^opt converted to the time domain using inverse FFT operations to produce FBE coefficients d^p^ in the time 

domain (step 412). The FBE coefficients b^pf are then applied with the FFE coefficients Qopt^o the DFE (as previously 
described with reference to FIG. 2. 
55 [0049] FIG. 5 IS a block diagram illustrating a discrete time fractionally-spaced DFE that operates according to the 
present invention. FIG. 6 is a block diagram illustrating a multhchannel equivalent of the discrete time fractionally- 
spaced DFE of FIG. 5. The approach used for the symbol-spaced model can be easily extended to fractionally spaced 
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models. In this section, we shall derive a fast algorithm for a 772-spaced equalizer. The fast solution for the general 
case of faster over sampling factors, say 773, 774, etc. will then follow simply by inspection of the arguments in this 
section. 

[0050] FlGs. 5 and 6 illustrate two equivalent representations for a discrete time model of a 772-spaced equalizer. 

5 FIG. 5 illustrates a DFE 500 in terms of a down sampler and an up sampler device, while FIG. 6 illustrates a DFE 600 
is a corresponding multichannel representation. Thus, as contrasted to FIG. 1 , the DFE 500 of FIG. 5 includes a down 
sampler 502 and an up sampler 504 that alter the sampling frequency. Further, the DFE 600 of FIG. 6 includes a first 
channel function 602 having channel estimate Hq and con-esponding noise function Vq (n) and FFE 606 coefficients Qq. 
Likewise, the DFE 600 of FIG. 6 also includes a second channel function 604 having channel estimate and corre- 

10 spending noise function (n) and FFE 608 coefficients [Qq, g^), 

[0051] The equivalence between the structures of FIG. 5 and FIG. 6 can be easily verified by writing each of the 
quantities {h,g] and v(n) in temns of their polyphase components, and interchanging their position with the down sampler 
and the up sampler (see, e.g., P.P. Vaidyanathan, "Multirate Systems and Filter Banks." Prentice Hail, NJ, 1993 and 
J.R. Treichler, I. Fijalkow, C.R. Johnson, Jr., "Fractionally spaced equalizers," IEEE Signal Processing Magazine, vol, 

15 13, no. 3, May 1996 for details). The quantities {h^ , h^] and [Qq, g^] are the so-called polyphase components of the 
channel and FFEs, and are denoted by A/ and L size vectors {Hq, h^,] and {gQ, g^} (see FIG. 6). The noise sequences 
{Vq (n), v^{n)) are also the even and odd samples of v(n), with powers 2a^. 

[0052] A model for the multichannel representation, just like the one forKhe symbol spaced case in equation (4) may 
be rewritten by collecting {go.Pil j^to a single vector g\ as: 

20 

flf'^tPi fl'ol Equation (44) 

with the output of the multichannel system given by: 

25 

[/o.n yi.nlflf' Equation (45) 

Thus, the following model of Equation (46) works for the input to g\ 

30 

[yo.n ^i.n 1 = 1^0 «il + l\n ^i.nl Equation (46) 

where [Hq, H^] are the convolution matrices associated with the subchannels {Hq, h^]. It is more convenient, however, 
^ to express the inner product in Equation (45) as a function of the original FFE coefficients vector g, by reordering the 
entries of g\ In this case, Equation (46) is replaced by: 

y„ = x^H'+ v„ Equation (47) 

40 

where 



55 



h(\) 



m 

h{2) 



h{\) 



h{2) 



h(2N-\) h{2N'2) 



KO) 
h(2) 

h{2N-\) h(2N-2) 



Equation (48) 



Given this convolution matrix and the noise variance a , the solution for the fractionally space problem is simply given 
by Equations (26) and (27), with {Hg ,H,o^} } replacedljy [H , H\ a^} (similariy to the symbol-spaced case, here we 
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assume M > A/ - 1). 

[0053] For the symbol-spaced DFE, the shift structure of implies that = ki^.^ „.^ (1 : L), That is, the normalized 
gain k[^ „ can be computed fast, by performing the order update of /c/^ to /ft+Ln-ii and then retaining the first L entries 
of kf^^^ as k^„. The procedure to derive a fast recursion In the T/2-spaced case follows the same principle. The only 
difference here is that the data matrix has now a double shift structure, in the sense that each row is formed by 
shifting two samples of the channel at a time. Thus, two successive order updates from /c^ ^.^ to /fi.+2,n-i are performed 
and then the first L entries of ^u2.im are retained: 



10 



Equation (49) 



so that 



15 



Equation (50) 



In other words, the resulting algorithm is given by two consecutive forward prediction problems similar to operations 
(1 ) through (6) of Table 2 of orders L and L+^, respectively. Table 3 lists the resulting algorithm for the 772-spaced case. 
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Initialization 



~ ^+1.-1 ^ ^L,0 ^ " ^ 

rt(0)=yi.,(0)=i 



For n=OtoS repeat operations (Hi to (\£S: 
(DForq^LtoL+l repeat operations (l)to(6\: 

(1) a^in-l) = hi2n + q-L+l)-h,.,wi^., 

(2) /,(/.-l) = ;',(«-l)a,(«-l) 

(4) - 1) = C((n - 2)+al(n - l)/,(n -1) 

(5) wi.-, =*<.,.j +k,^.J,(n-l) 

(n)*,^=W.a:^) 



Table 3: Fast Transversal Computation of FFE coeflidents for the r/2-spaced equalizer. 



^ [0054] In the general case of a 77S spaced equalizer, the fast DFE tap computation is a straightfonvard extension 
of the above algorithm, where q is iterated from L to (L+S-1). However, in the fractionally spaced equalizer, the RLS 
problem is formulated as a multi-channel problem. In such case, a multi-channel Kalman gain is calculated for the 
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multi-channel RLS problenn and the FFE taps are deternnined there from. Note that in the expressions of Table 3, 

successive order updates are perfonned. 

[0055] Now, the optimal FBE taps could now be computed as 



where F^^F® Ig and Is a block circulant matrix, which is formed by extending the matrix H in the same manner 
as in Equation (39). As was the case with the symbol spaced operations described above, the convolution may be 
perfonned in a transfomied domain, e.g.. Frequency Transfomnation domain, Discrete Cosine Transformation domain 
and the Discrete Hadamard Transformation domain, among others. In such case a Croniker product may be employed 
in conjunction with such domain transfonnatlon in perfomning multi-channel convolution In the selected transformed 
domain. 

[0056] FIG. 7 is a block diagram illustrating a transceiver constmcted according to the present invention. The com- 
ponents of the transceiver 700 are resident in a communications device and are illustrated only generally to show how 
the operations of the present invention would be accomplished in such a transceiver. The transceiver 700 includes a 
receiver section 702 and a transmitter section 704 and further includes either a wireless transceiver 706 or a wired 
transceiver, depending upon the system in which the transceiver 700 is implemented. In the case of a cellular device, 
RF device, satellite system device, or another wireless implementation, the transceiver 700 includes a wireless trans- 
ceiver. However, in the case of a cable modem, a LAN device, a home networking device, or another device that 
couples to a physical media, the transceiver 700 includes a wired transceiver 708. Further, if the present invention is 
implemented in a Seriatizer/Deseriaiizer (SERDES) or similar application, the receiver section 702 and transmitter 
section 704 may couple to a media without a wired transceiver 708. 

[0057] Further discussion of the transmitter section 704 and receiver section 702 are in the context of baseband 
processing. In such case, the receiver section 702 receives a baseband signal from the wireless transceiver 706 (or 
wired transceiver 708) that is baseband modulated and operates upon the baseband signal to extract data. These 
operations include detenmining DFE coefficients according to the present invention and operating upon the baseband 
signal using the determined DFE coefficients. 

[0058] The transmitter section 704 receives digital data to be transmitted from a host, codes the digital data into a 
baseband signal, and passes the baseband signal to the RF transceiver 706. The RF transceiver 706 couples the 
baseband signal to an RF carrier to create an RF signal and transmits the RF signal to a receiving device across a 
wireless link. 

[0059] The receiver section 702 receives a baseband signal that carries coded data from the RF transceiver 706. A 
Programmable Gain Amplifier (PGA) 712 adjusts the gain of the baseband signal and then provides the gain-adjusted 
baseband signal to an Analog-to-Digital Converter (ADC) 714 for sampling. The ADC 208 samples the gain adjusted 
baseband signal at a particular sampling frequency, (that is the symbol clock frequency), to produce samples thereof. 
[0060] A processor 71 0 couples to the output of the ADC 71 4 and analyzes a preamble sequence contained in each 
received physical layer frame. Based upon the preamble sequence, the processor 710 detenmines a gain to be applied 
to portions of the baseband signal corresponding to the data carrying portions of the frame and provides this gain to 
the PGA 71 2. Further, the processor 71 0 may also interact with the optional timing compensation section 71 6 to com- 
pensate for symbol timing and RF carrier mismatches. 

[0061] The processor 710. based upon the preamble sequence (and based upon actual extracted data in some 
operations), also determines FFE 104 and FBE 108 coefficients. The manner in which these coefficients are detemiined 
was previously described in detail herein. Further, and as was also previously described, the processor 710 may es- 
timate a channel and calculate DFE coefficients based upon unknown but assumed data content. After the processor 
71 0 detemnines these coefficients, they are applied to the FFE 1 04 and FBE 1 08 for subsequent use in extracting data 
from the baseband signal. 

[0062] The stnjcture described in FIG. 7 may be embodied using various types of circuits formed using various 
manufacturing processes. For example, in one particular embodiment, the RF transceiver 706 (or wired transceiver 
708) is embodied in a first integrated circuit that is coupled to a second integrated circuit that includes the transmitter 




Equation (51) 



where is a block diagonal matrix, which satisfies 



Equation (52) 
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section 704 and the receiver section 702, among other circuits. In another embodiment, the RF transceiver 706, the 
transmitter section 704 and the receiver section 702 are all formed on a single monolithic integrated circuit. These 
integrated circuits may be constructed in CMOS or another semiconductor technology, e.g., PMOS, NMOS, Bipolar, etc. 
[0063] Further, the receiver section 702 of FIG. 7 may be constructed using various circuit elements/combinations. 

5 In one embodiment, all structures past the ADC 714 in the receiver section 702 are embodied using a Digital Signal 
Processor (DSP) or similar processing device. In another embodiment, dedicated signal path circuitry embodies each 
of the structural components of the receiver section 702, including the processor 710. While a DSP implementation 
would provide more flexibility, dedicated signal path circuitry would typically provide higher performance at a lower cost 
and with lower power consumption. 

10 [0064] The structure and operation of the present invention may be employed in satellite comnriunication systems, 
satellite television systems, HDTV systems, fixed wireless communication systems, mobile wireless communication 
systems, cable modem/television communication systenns, home networking systems, wireless local area networking 
systems, wired local area networking systems, and many other types of communication systems. The present invention 
applies to all types of communication devices in which equalizers are employed to operate upon received signals. 

15 [0065] FIG. 8 is a block diagram illustrating a Multi-Input-Multi-Output (MIMO) digital communication system that 
operates according to the present invention to equalize a channel. MIMO digital communication systems, as the temi 
Implies, are communication systems that include multiple inputs on a transmitting side and multiple outputs on a re- 
ceiving side. In such systems. MIMO decision feedback equalization is used to mitigate inter-symbol Interference (ISI) 
that results from channel multi path propagation. 

20 [0066] In the embodiment of FIG. 8, an input symbol stream includes P unique transmitted signals, represented by 
XQ(n), x^{n), .... Xp.i(n) (as are further shown with reference to FIGS. 8, 9, and 10). The nomenclature of FIG. 8 is used 
herein to describe the method of the present invention. Thus, the P transmitted signals are referred to in combination 
as transmitted signal vector x{n). The transmitted signal vector x(n) consists of a known training sequence followed by 
unknown data. The transmKted signal vector x(n) passes through a channel 802 represented by H{z) that is represented 

25 to have N taps. The channel 802 includes additive noise v(n), which Is white and Gaussian and has power spectral 
density a^. The output of the channel 804 is referred to as y(n) and has a length O (the value Q is a function of P and N). 
[0067] \ MIMO DFE that equalizes the channel includes a MIMO FFE 804 having function G(z) and length L and a 
MIMO FBE 808 having function B(z) and length M (that is assumed to be greater or equal to the channel memory, i. 
e., M > A/-1). The outputs of both the MIMO FFE 804 and the MIMO FBE 808 have a width P that corresponds to the 

30 multiple numbers of signals. MIMO decision block 806 makes soft decisions for the MIMO symbol stream and produces 
as P output symbol streams jf o(n - 6). x iC'^ - 5)i x p.^('^• S)- These soft decisions are provided to further receiver 
components, e.g., a Viterbi decoder to produce hard decisions based upon the soft decisions. 
[0068] Examples of systems in which the MIMO DFE of the present invention may be implemented include wireless 
systems, e.g., cellular wireless systems, fixed local loop wireless systems, wireless local area networks, high definition 

35 television systems, etc., and wired systems in which multiple transmitters and receives share a common media, e.g., 
cable modem systems, etc. Diagrams representing such systems are illustrated in FIGs. 9, 10, and 11. The structure 
of FIG. 7 may be employed to implement the structure and operations described with reference to FIGs. 8-11. 
[0069] FIG. 9 is a system diagram illustrating a wireless digital communication system 900 in which a MIMO receiver 
904 operates according to the present invention. The wireless system 900 includes a MIMO transmitter 902 that op- 

^0 erates to receive a plurality of signal streams and to form a wireless signal and that transmits the wireless signal. The 
MIMO receiver 904 receives the wireless signal after a corresponding channel has operated upon it. The MIMO receiver 
904 includes a MIMO DFE that operates according to the present invention to equalize the received signal and to 
produce output based thereupon. As is illustrated, the MIMO receiver 904 may include antenna diversity structure. 
[0070] FIG. 1 0 is a system diagram illustrating a wireless digital communication system 1 000 that includes a plurality 

45 of transmitters 1 002A-1 002G and a MIMO receiver 1 004 that operates according to the present invention. Each of the 
plurality of transmitters 1002A-1002G produces a respective transmitted signal that carries a respective input symbol 
stream, x,(n), X2(n), .... Xp(n). The MIMO receiver 1004 receives the plurality of transmitted signals as a composite 
signal after they have been operated upon the channel, equalizes the plurality of received signals, and produces soft 
decisions for the equalized symbols. The wireless digital communication system 1000 of FIG. 10 may be embodied as 

50 a cellular wireless network, a fixed wireless access network, a wireless local area network, or in another type of wireless 
communication system in which a MIMO receiver 1004 receives wireless transmissions from a plurality of wireless 
transmitters 1002A-1002G. As Is illustrated, the MIMO receiver 904 may include antenna diversity structure. 
[0071] FIG. 11 is a system diagram illustrating a wired digital communication system 1100 that includes a plurality 
of transmitters 1 1 02A-1 1 02G and a MIMO receiver 11 04 that operates according to the present invention. Each of the 

55 plurality of transmitters 11 02A-1 1 02G produces a respective transmitted signal that carries a respective input symbol 
stream, XQ(n), Xi(n), Xp.,(n) and couples its respective signal to a network infrastructure 1106. The MIMO receiver 
1 1 04 couples to the network infrastructure 1 1 06 and receives the plurality of transmitted signals as a composite signal 
after they have been operated upon the channel, equalizes the plurality of received signals, and produces soft decisions 
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for Ihe equalized symbols. The wired digital communication system 1 000 of FIG. 1 1 may be embodied in a cable modem 
system, or In another type of system in which the MIMO receiver 1004 receives wired transmissions from a plurality 
of transmitters 1002A-1002G 

[0072] Referring again to FIG. 8, an efficient equalization technique for equalizing a channel in a MIMO system 
includes first estimating the channel impulse responses between each transmitter and each receiver using the training 
sequence, and then using this estimate to compute the optimal DFE tap coefficients corresponding to the estimated 
channel. The computed tap coefficients are then uploaded to the equalizer taps of the MIMO FFE (including MIMO 
FFE 804 and MIMO FBE 808) of FIG. 8. This procedure should be repeated as often as possible, especially in cases 
of fast varying channels. Moreover, the received data stream is usually buffered during the period assigned for the 
channel estimation and the equalizer tap computation. 

[0073] In this context, a critical factor for the success of this equalization structure is the complexity of the MIMO 
DFE tap computation. A fast computation technique has the following benefits: 

1 . It reduces the memory size required to buffer the received sequence during the period required for tap compu- 
tations. 

2. It allows more frequent uploading of new equalizer coefficients, thus enabling the equalizer to track fast channel 
variations. 

3. It simplifies the needed hardware, especially If such computations are perfonned through structured recursive 
equations. 

[0074] According to the present invention, an efficient method for computing the optimal MIMO MMSE-DFE coeffi- 
cients includes characterizing a solution to the problem of solving for MIMO FBE and FFE coefficients as a fast recur- 
sive-least-squares (RLS) adaptive algorithm. This fast algorithm has the following advantages over the current most 
efficient algorithm: (1 ) the proposed algorithm has less overall computational complexity than prior techniques in solving 
for DFEs having the same channel and DFE filters length; and (2) the technique relies on the use of a set of structured 
recursions, which makes it attractive in data-path implementations which do not rely on the use of a DSP core. 
[0075] For simplicity of presentation a symbol-spaced MIMO DFE is considered. Extension of the teachings provided 
herein to the fractionally-spaced case is straightfonvard. Further, the following commonly used assumptions are made 
herein: 

1 . The input sequences to the /-th channel, {x, (n)}, are complex, independent and identically distributed with power 
. The autocorrelation matrix of the corresponding vector process x{n) is given by: 



R« = /<8><//a9(aJ,o*^MM..,<p.,) Equation (53) 

where ® denotes the Kronecker product, and I is identity matrix of appropriate dimensions. 

[0076] The noise sequences [v^ (n)] are assumed to be white Gaussian with power o^^ . The autocorrelation matrix 
of size LQ of the corresponding vector noise process v(n) is given by: 



Rv = ® diag (a,%, a^,, a^Q,,) Equation (54) 

A 

[0077] The decisions (from decision block 806) x(n - 5) are assumed to be correct, and hence equal to x(n - 5) . 
Moreover, the assumption that only previous decisions of other users are available at the present time, I.e.. time (n - 
5). It is further assumed that the use of current decisions of other users contributes with negligible Improvement of the 
final output SNR. 

[0078] The MIMO FFE 804 has a matrix filter G(z) with length L. The number of matrix taps Mof the MIMO FBE 808 
having matrix filter B(z) Is such that M^N-^. 

[0079] The goal in detenmining the taps of the MIMO FFE 804 and the MIMO FBE 808 is to minimize the mean 
squared error quantity: 

^ = E||x(n - 5) - x(n - 5)f , Equation (55) 
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where x(n - 6) is the delayed input signal estimate prior to the decision. By collecting the tap nnatrix coefficients both 
G(z) and B(z) into vectors, the received signal x(n - 5) is expressed as: 



x(n - 5) = y„g - x„b, 



Equation (56) 



where: 



10 



Equation (57) 



[0080] H is the (A/ + L - 1 )P X LO convolution matrix associated with the channel matrix coefficients h(0), h(1 ) h 

(A/ - 1), and is given by: 



15 



H: 



25 



h(0) 

h(l) 



0 
h(0) 



hiN-l) h(N-2) 
0 h(A^-l) 



0 
0 

h(0) 
h(l) 

h(Ar-l) 



Equation (S8) 



The row vector x^ is the 1 x(N+ L-1)P input regressor: 



30 



Equation (59) 



With x(n) A [xo (n) (n) ... Xp,^ (n)]. The row vector y„ is the 1 x LQ Input regressor to the feed fonward filter g. I.e., 



35 



yAbin) y(n-l) y(«-i+l)] 



Equation (60) 



40 



Equation (61) 



45 



Where y(n) ^ l Vb (") Vi i^) Xq-i ("P and (gj) has size Q x P. Similarly, x „ is the 1 x MP input regressor to the 
(strictly causal) feedback filter b, i.e., 



Equation (62) 



50 



col[bo,b,,...,b;^^.J 



Equation (63) 



Where [bj are P x P. Also, v^ Is the 1 x LQ noise vector process. 
55 [0081] By collecting g and b into a single matrix w, the minimization of Equation (55) can be written as: 
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minE 



x(/i-J)-[y„-xJ 



Equation (64) 



[0082] Now, denoting as the variance of the augmented Input regression vector u, and Rux(n-6) as the cross var- 
iance of the augmented input regression vector u, a solution to this smoothing problem may be determined using a 
10 |<nown technique, see, e.g., T. Kailath, A.H. Sayed, and B. Hassibi, Linear Estimation, Prentice Half, N.J., 2000. to 
produce an optimal solution for w according to: 



15 



where: 



Wopt - ^ux(n-6)» 



Equation (65) 
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Equation (66) 



and 



25 



30 



Equation (67) 



[0083] Using a channel output model, such as the channel output model described in T. Kailath, A.H, Sayed, and B. 
Hassibi, Linear Estimation, Prentice Hall, N.J., 2000, and the fact that {x,{n)} are individually Identically distributed (i.i. 
d.), The following closed form expressions for 



35 



may be formed: 



40 



Ry=Ey„y„ =R^ +H*R^H. 



Equation (68) 



45 



R^=H(ExX) = H- 



Equation (69) 



50 



Rj — Rj 



Equation (70) 



R,.(„-#)-H*Ex;x(n-^) = H- 



o 



Equation (71) 
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Equation (72) 



where H is a submatrix of H, such that 



H = 



10 



H 

1^2J 



, where H, is (S + l)PxLQ 



Equation (73) 



[0084] Now, with the quantities defined above, Equation (65) becomes: 



15 



'[ -S.H R. J L O 



Equation (74) 



20 Using the well-known inverse of block matrices, Equation (74) may be rewritten as: 



25 



30 



= [I] (R, + KR<"H. + H;Ri,=»H,)-'h'R';' . 
[0085] Further, because M ^ A/ -1 , the quantities hg, and H2 are such that: 



Equation (75) 



Equation (76) 



This implies the following expressions for the optimal feedback and feed forward coefficients may be rewritten as: 



35 



40 



g,p,=(R,+H:R</JH,)»h;R<;> 



Equation (77) 



Equation (78) 



where R^^^ and R^'^ are size 6 and P matrices. The above expressions are valid for all values of the delay 5. In general, 
45 the optimiil value for the decision delay 5 Is within the range 



L-1<6opt^A/+l.-2. 



50 [0086] In finding a solution for g^p,, a coefficient matrix is first determined as: 



Equation (79) 



55 



Where the definition Hg^ R^^Hg is made. The optimal solution for the feed fon/vard coefficients is then given by: 
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9opt = P8.Di9i'^- Equation (80) 

[0087] A reader acquainted with the theory of recursive (east squares algorithms might readily recognize that the 
quantity gopt = ^ail^ conresponds exactly the definition of the Kalman gain matrix in the recursive least squares algo- 
rithm. The Kalman gain matrix is used to update the optimal weights in a regularized (multi channel) RLS problem. 
More specifically, given a (n + 1)P x LO data matrix H„ and the corresponding coefficient matrix the Kalman gain 
g„ = P^ h can be time-updated according to the following recursions (see e.g. T. Kaliath, A.H. Sayed, and B. Hassibi, 
Linear Eslimation, Prentice Hall, N.J., 2000): 

y'^ (n) = Ip +h^P„.^ ,h^, Equation (81) 

9n = Pn-1 t!„Y('')' Equation (82) 

Pn =Pn.r9nY'^ (n)gn. Equation (83) 

where P.1, =R'^ and go = 0. 

[0088] The computation of g „ can be done efficiently, via a fast RLS recursion. That is, consider Equation (82), which 
can be written as: 

Qn = M(n), Equation (84) 

where y(n) is defined in Equation (81). The quantity k„ = P^ih , is referred to as the norma//zed Kalman gain matrix 
in adaptive RLS filtering. The key to achieving fast recursions is {o propagate k„ efficiently, by successive order updates 
and order down dates of k„, using forward and backward LS prediction problems. The well-known fast transversal filter 
(fH'Fj is an example where such fast technique Is encountered, see e.g., J. Cioffi and T Kailath, "Fast recursive-least- 
squares transversal filters for adaptive filtering," IEEE Trans, on Acoust.. Speech Signal Processing, vol. ASSP-32, 
pp. 304-337, April1984. 

[0089] The difference here, however, is that the entries of the regressor are the channel matrix coefficients h(0),..., 
h(A/- 1), each of size P x O (In the multi channel RLS equations, the entries of the regressor are usually given by a 
row or a column vector, and the recursions have no matrix inversions). 

[0090] Now let Kq „ = k„, and denote {k/^^} the augmented Kalman gains that correspond to estimating each column 
of h(n). say. [h(n)], / from [h(/?)] ;,/4.io!!n-i- order to update k^ „ efficiently without matrix inversions, we perform Q 
successive order updates from ko^i to kQ^^. Now, note that in general the fast recursions that computer kQ rely on 
fonvard and backward prediction problems. Here, because M > N-1 the matrix Hg has a block lower triangular structure, 
the desired signal for the backward prediction problem is always equal to zero. As a result, the least-squares backward 
prediction solution up to time 6 will be equal to zero, and the fast algorithm will need only the forward prediction part. 
That is, 

'fo.n-i „.i ->...->kQ„.i Equation (85) 

so that 

ko.n =ko.n.i • Equation (86) 

[0091] The resulting algorithm is then given by consecutive updates of 0 forward prediction problems where we 
estimate the columns of h(n) from h^, and compute the corresponding forward prediction matrix-valued quantities. 
Table 4 lists the resulting multi channel algorithm. 
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Initialization 

For a = 0(0 0-1 perform the following initialization steps: 

ro(0)=i. 



For n = OtoS. repeat operations flP. and (ISi : 
(I) Forq = 1 to Q. repeat CH to ffl : 

(1) a,(«-l) = [h(«)lo^-h^,w^,., 

(2) f,(«-l) = y,(Af-l)a,(«-l) 

(4) («-!) = C{{n - 2) + o; (« - l)f,(n - 1) 

(5) w^,..=<^.,+k,..^.,f,(«-l) 

(ID k,,=ke„.,(l:ei.:) 
(ni)yo(«) = re(»-l) 



Finally, set g^=k,^y,(.5)R'.' 



1/2 



Table 4: Fast Computation of for the MIMO DFE. 
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[0092] The exponentially weighted FTP algorithm can become unstable in some finite precision implementations. 
Hence, the reader may wonder whether this is also the case for the present simplified algorithm, since both have the 
same essence. There aire two advantages, however, that prevent the simplified algorithm from becoming unstable: 

(i) The main source of error propagation in the full FTP algorithm arises in the backward prediction section of its 
recursions. Here, by ruling out the equations associated with the backward prediction problem many of the recursive 
loops that contribute to the unstable behavior of the full FTP algorithm are automatically eliminated. 

(ii) Another source of instability of the FTP recursions is related to the forgetting factor X that appears in the expo- 
. nentially weighted RLS problem. Theoretically with X < 1 , the redundant components generated by numerical errors 

would decay to zero as A/-> «>. However, an averaging analysis shows that this will lead to unstable modes at 1/X 
that will cause instability in finite precision. Now, note that the present fast recursions deal with the problem of 
filtering a finite set of data samples of the channel model. In other words, for the purposes of the present invention, 
the algorithm must stop when n = 6. Moreover, in the present application, the con-esponding forgetting factor is 
always equal to one, in which case the recursions will present much better numerical behavior. This means that 
even if the simplified fast algorithm were to become unstable, it would not be likely to happen within the first 5 
iterations. Purthemrv^re, if this were still the case, a simple increase of the word length would overcome the problem. 

[0093] The computation of b^pj in equation (36) can be performed efficiently via using fast Fourier transform tech- 
niques. This can be done by extending the block Toeplitz structure of R to form aKPxKQ block-circulant matrix CpQ, 
where each block entry has size P x 0(see, e.g., M. Vollmer, J. Gotze, M Haardt. "Efficient Joint Detection Techniques 
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for TD-CDMA In the Frequency Domain," COST 262 Workshop "Multiuser Detection in Spread Spectrum Communi- 
cations", Ulm, Germany, 2001 ). 

[0094] In the present case, there is a block diagonal matrix ApQ such that CpQ = F^Ap qFq where F/ are block- 
Fourier transforms given by F; = FS>/j,with 0 denoting the Kronecker product and the matrix Fthe K x KDFT matrix. 
The diagonal elements of ApQ can thus be computed from the block-DFT of the first block-column of OpQ. The bopi 
then becomes 



Hence, the complexity is given b P inverse FFTs, O FFTs, and P FFTs of size K, and KPQ complex multiplies. Thus 
the overall complexity for computing the feedback filter is KPQ + (2P + 0)Klog2K). 

[0095] FIG. 12 is logic diagram illustrating operation according to another embodiment of the present invention for 
determining FFE coefficients. The operations of FIG. 1 2 differ from the operations of FIG. 3 in that the delay of FBE is 
not limited to being less than or equal to the length of the FFE but is chosen to correspond to the channel being 
equalized (step 1 202). FIGs. 1 4 and 1 5 describe differing techniques for selecting the DFE delay. Referring still to FIG. 
1 2, operation continues by formulating a solution that, when solved, will yield the DFE coefficients, wherein the solution 
is formulated as a least squares problem that is based upon the channel response (step 1204). 
[0096] Next, the operation considers whether the DFE delay is greater than the length of the FFE minus one, i.e., is 
the DFE delay greater than the FFE delay (step 1206). When the delay of the DFE is less than or equal to a number 
of taps of the FFE minus one, the method includes solving the least squares problem to yield FFE coefficients of the 
DFE coefficients by performing a first set of operations (step 1 208). However, when the delay of the DFE is greater 
than the delay of the FFE, operation includes solving the least squares problem to yield FFE coefficients of the DFE 
coefficients by performing a first set of operations (step 1210) and by performing a second set of operations (step 
1212). Finally, the method includes (from step 1208 and from step 1212) convolving the FFE coefficients with a con- 
volution matrix that is based upon the channel response to yield FBE coefficients of the DFE coefficients (step 1214). 
[0097] The operations of FIG. 12 are further described in the following and Illustrated in Table 5. Tlie solutions pro- 
vided by Tables 2. 3 and 4 were valid only for 8 ^ L - 1 , which Is inadequate for channels with long pre-echoes. The 
general solution for the FFE taps was indicated in Table 1 . However, for the values of n < L - 1 , the variables corre- 
sponding to the backward prediction problem are ail equal to zero, while the fonward prediction variables are not. For 
the values of n > L - 7, the backward prediction variables are non-zero as well. To exploit this property to reduce 
complexity of the full FTF algorithm in Table 1 , updating the backward prediction variables forn > L - 7 is done, as follows: 




MP 




Equation (87) 
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Initialization 

7'«(0) = 1 

For n^OtoL'l repeat operations (H to (7): 
(2)f(n-l) = r,in-\)a{n-l) 

(4) (« - 1) = ^/ (n - 2) + a * (« - l)/(n - 1) 

(5) h/,=m/,+A,,,V(«-1) 

0) ^L.n ~K.L.n-X 
For n = Lto S repeat operalions CI) to fl3): 



a,{n-\) = h{n)-h„_,wi.^ 
fin-\) = r,{n-\)a{n-\) 

C^(B-l) = ^^(«-2) + a*(«-l)/(«-l) 
1*^-1 = '^-2+*i..-i/(n-l) 

v(n) = (last entry of A^^, ,) 

=*.:i.»+v(«)v»^., 
A«) = ^»(/i-l)v*(«) 

" i-f...(«)A(«K(«) 
K«)=rt(«)A(«) 

C*(n) = C*(«-l) + /9*(«)6(«) 



(1) 

(2) 
(3) 

(4) 
(5) 

(6) 

(7) 
(8) 
(9) 

(10) 

(11) 
(12) 
(13) 

If n = «y set g^^kgY^aS) 



Table 5: FTF computation of FFE taps for long equalizer delays 
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[0098] By so doing, in the first L iterations, we need 3172 multiplies/iteration. For the last (5-L) Iterations, the complexity 
is 5172 multipltesAiteration. Thus the total complexity is [3L2 + 5L (S-L)]/2 multiplies. If the full FTF algorithm was used, 
we would need 5L6/2, which is significantly more complex. 

[0099] FIG. 1 3 is a graph illustrating a hypothetical channel impulse response indicating how channel delay may be 
determined according to the present invention. To determine the value of the delay 5 of the FBE, the following procedure 
could be used. The channel delay 5^ could be computed as the index of the strongest channel tap. The equalizer delay 
could then be computed as 5=5c+L-1 . 

[0100] FIG. 14 is a logic diagram illustrating operation according to a first embodiment of the present invention in 
selecting FBE coefficients. This operation includes first selecting an initial equalizer delay value (step 1402). Then, for 
the particular equalizer delay value (FBE delay), the operation includes calculating a set of FFE coefficients and FBE 
coefficients (step 1404). Then operation then includes determining the FBE tap energy for the particular FBE delay 
(step 1406). If this FBE tap energy is less than a FBE tap energy threshold (as determined at step 1408) the FFE and 
FBE coefficients determined for the particular FBE are employed (step 1414). If the FBE tap energy is not less than 
the FBE tap energy threshold (as detennined at step 1 408), it is next determined whether the FBE delay is at a maximum 
(step 1 41 0). If so. no acceptable FBE coefficients have been determined for the particular channel response and FBE 
configuration, the operation has failed (step 1416). However, if the FBE delay is not at a maximum (as detemnined at 
step 1410), the equalizer delay is increased by a step size (or othenvise altered as may be appropriate for the given 
operations) at step 1412 and operation returns to step 1404. 

[0101] FIG. 15 is a logic diagram illustrating operation according to a second embodiment of the present invention 
in selecting FBE coefficients. This operation includes first selecting an initial equalizer delay value (step 1502). Then, 
for the particular equalizer delay value (FBE delay), the operation Includes calculating a set of FFE coefficients and 
FBE coefficients (step 1504). Then operation then Includes determining the FBE tap energy for the particular FBE 
delay (step 1506). The operation then Includes storing the FFE coefficients, the FBE coefficients, and the FBE tap 
energy along with the particular FBE delay (step 1508). It is next detemiined whether the FBE delay is at a maximum 
(step 1510). If so, the FFE coefficients, FBE coefficients, and FBE delay is selected that corresponds to a minimum 
FBE tap energy (step 1516) of all solutions that have been stored. If the FBE delay is not at the maximum FBE delay 
(as determined at step 1510), the equalizer delay is increased by a step size (or otherwise altered as may be appropriate 
for the given operations) at step 1512 and operation returns to step 1504. 

[0102] The invention disclosed herein is susceptible to various modifications and alternative fomns. Specific embod- 
iments therefore have been shown by way of example in the drawings and detailed description. It should be understood, 
however, that the drawings and detailed description thereto are not intended to limit the invention to the particular form 
disclosed, but on the contrary, the Invention Is to cover all modifications, equivalents and alternatives falling within the 
spirit and scope of the present invention as defined by the claims. 



Claims 

1. A method for computing Decision Feedback Equalizer (DFE) coefficients, the method comprising: 

estimating the channel response of a channel operated upon by the DFE; 

formulating a solution that, when solved, will yield the DFE coefficients, wherein the solution is formulated as 
a least squares problem that is based upon the channel response; 

solving the least squares problem to yield Feed FonA/ard Equalizer (FFE) coefficients of the DFE coefficients, 
wherein solving the least squares problems Includes: 

selecting a delay of the DFE; 

when the delay is less than or equal to a number of taps of the FFE minus one, performing a first set of 

operations; and 

when the delay Is greater than the number of taps of the FFE minus one, perfonming a second set of 
operations; and 

convolving the FFE coefficients with a convolution matrix that is based upon the channel response to yield 
Feed Back Equalizer (FBE) coefficients of the DFE coefficients. 

2. The method of claim 1 . wherein the least squares problem is formulated as a Kalman gain solution. 

3. The method of claim 2, wherein the Kalman gain solution is detenmined using a Fast Transversal Filter (FTF) 
algorithm. 
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4. The method of claim 3, wherein when the DFE delay is chosen to be less than or equal to a number of taps of the 
FFE minus one, a length of the F6E is chosen to force the FTF algorithm to use a lower triangular matrix. 

5. The method of claim 1 . wherein determining the delay of the DFE includes: 

5 

determining a delay of the channel based upon the channel response of the channel; and 

adding the delay of the channel to the number of taps of the FFE minus one to produce the delay of the DFE. 

6. The method of claim 5, wherein the delay of the channel Is selected as the length of the precursor of the channel 
10 response. 

7. The method of claim 1 . wherein determining the delay of the DFE includes: 

calculating a plurality of sets of FFE coefficients and FBE coefficients, wherein each of the plurality of sets of 
15 FFE coefficients and FBE coefficients corresponds to a respective DFE delay; and 

selecting one of the plurality of sets of FFE coefficients and FBE coefficients based upon the tap energy of 
the respective FBE coefficients. 

8. The method of claim 7, wherein selecting one of the plurality of sets of FFE coefficients and FBE coefficients based 
20 upon the tap energy of the respective FBE coefficients comprises selecting a set of FFE coefficients and FBE 

coefficients that meets an FBE coefficient tap energy threshold. 

9. The method of claim 7, wherein selecting one of the plurality of sets of FFE coefficients and FBE coefficients based 
upon the tap energy of the respective FBE coefficients comprises selecting a set of FFE coefficients and FBE 

25 coefficients tiiat provides a minimum FBE coefficient tap energy. 

10. The method of claim 1, wherein in convolving the FFE coefficients with the convolution matrix that is based upon 
tfie channel response to yield the FBE coefficients of the DFE coefficients, the method further comprises: 

30 extending a convolution matrix created based upon the channel response to a bigger circulant nr^atrix; and 

perfomning the convolution in a transformed domain. 

11. The method of claim 10, wherein the transformed domain is selected from the group consisting of the frequency 
domain, the Discrete Cosine Transform domain and the Discrete Hadamard Transform domain. 

35 

12. The metiiod of claim 1, wherein in convolving the FFE coefficients with tiie convolution matrix that is based upon 
the channel response to yield the FBE coefficients of the DFE coefficients, tiie method further comprises: 

extending a convolution matrix created based upon the channel response to a bigger circulant matrix; and 
^0 computing the convolution in the frequency domain. 

13. The method of claim 12, wherein computing tiie convolution in the frequency domain Includes: 

transfomning the convolution matrix and the FFE coefficients from the time domain to the frequency domain 
using a Fast Fourier Transfonm; 

computing the convolution in the frequency domain to produce the FBE coefficients; and 
transfomning the FBE coefficients from the frequency domain to the time domain. 

14. The method of claim 1 , wherein the channel response is based upon a known preamble sequence of a packet 
50 upon which the DFE operates. 

15. The metiiod of claim 1 , wherein the channel response is based upon decisions made for data carried in a packet. 

16. The metiiod of claim 1 , wherein the channel response is based upon: 

55 

a known training sequence contained in a packet preamble; and 
decisions made for data carried in the packet. 
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17. The method of claim 1 , wherein the FFE and the FBE of the DFE are symbol spaced. 

18. The method of claim 1 , wherein the FFE and the FBE of the DFE are fractionally symbol spaced. 

5 19. The method of claim 1 , wherein the method is applied within a Digital Television Broadcast system. 

20. The method of claim 1 , wherein the method is applied within a Satellite communication system. 

21. A method for computing Decision Feedback Equalizer (DFE) coefficients, the method comprising: 

10 

estimating the channel response of a channel operated upon by the DFE; 

formulating a solution that, when solved, will yield the DFE coefficients, wherein the solutton is formulated as 
a least squares problem that is based upon the channel response; 

solving the least squares problem using a Kalman gain solution to yield Feed Forward Equalizer (FFE) coef- 
15 ficients of the DFE coefficients using a Fast Transversal Filter (FTF) algorithm, wherein solving the least 

squares problenns includes: 

selecting a delay of the DFE; 

when the delay is less than or equal to a number of taps of the FFE minus one, perfonming a first set of 
20 operations; and 

when the delay is greater than the number of taps of the FFE minus one, performing a second set of 
operations; and 

convolving the FFE coefficients with a convolution matrix that is based upon the channel response to yield 
25 Feed Back Equalizer (FBE) coefficients of the DFE coefficients. 

22. The method of claim 21 , wherein when the DFE delay is chosen to be less than or equal to a number of taps of 
the FFE minus one, a length of the FBE is chosen to force the FTF algorithm to use a lower triangular matrix. 

30 23. The method of claim 21 , wherein detennining the delay of the DFE includes: 

determining a delay of the channel based upon the channel response of the channel; and 

adding the delay of the channel to the number of taps of the FFE minus one to produce the delay of the DFE. 

35 24. The method of claim 23. wherein the delay of the channel is selected as the length of the precursor of the channel 
response. 

25. The method of claim 21, wherein determining the delay of the DFE includes: 

40 calculating a plurality of sets of FFE coefficients and FBE coefficients, wherein each of the plurality of sets of 

FFE coefficients and FBE coefficients corresponds to a respective DFE delay; and 
selecting one of the plurality of sets of FFE coefficients and FBE coefficients based upon the tap energy of 
the respective FBE coefficients. 

45 26. The method of claim 25, wherein selecting one of the plurality of sets of FFE coefficients and FBE coefficients 
based upon the tap energy of the respective FBE coefficients comprises selecting a set of FFE coefficients and 
FBE coefficients that meets an FBE coefficient tap energy threshold. 

27. The method of claim 25, wherein selecting one of the plurality of sets of FFE coefficients and FBE coefficients 
50 based upon the tap energy of the respective FBE coefficients comprises selecting a set of FFE coefficients and 

FBE coefficients that provides a minimum FBE coefficient tap energy. 

28. The method of claim 21, wherein the FFE and the FBE of the DFE are symbol spaced. 

55 29. The method of claim 21 , wherein the FFE and the FBE of the DFE are fractionally symbol spaced. 

30. The method of claim 21, wherein the method is applied within a Digital Television Broadcast system, or within a 
Satellite communication system. 
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31. A method for computing Decision Feedback Equalizer (DFE) coefficients, the method comprising: 

estimating the channel response of a channel operated upon by the DFE; 

formulating a solution that, when solved, wilt yield the DFE coefficients, wherein the solution is formulated as 
a least squares problem that is based upon the channel response; 

for each of a plurality of FBE delays, solving the least squares problem to yield a set of Feed Forward Equalizer 
(FFE) coefficients and Feed Back Equalizer (FBE) coefficients of the DFE coefficients, wherein each set of 
FFE coefficients and FBE coefficients corresponds to a respective FBE delay; and 
selecting a set of FFE coefficients and FBE coefficients based upon the tap energy of the FBE coefficients of 
the set. 

32. The method of claim 31 , wherein selecting a set of FFE coefficients and FBE coefficients based upon the tap 
energy of the FBE coefficients of the set comprises selecting a set of FFE coefficients and FBE coefficients that 
meets an FBE coefficient tap energy threshold. 

33. The method of claim 31 , wherein selecting a set of FFE coefficients and FBE coefficients based upon the tap 
energy of the FBE coefficients of the set comprises selecting a set of FFE coefficients and FBE coefficients that 
provides a minimum FBE coefficient tap energy. 
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